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Abstract 

We comment on the new trend in mathematical physics that con- 
sists of obtaining Taylor series for fabricated linear and nonlinear un- 
physical models by means of homotopy perturbation method (HPM), 
homotopy analysis method (HAM) and Adomian decomposition method 
(ADM). As an illustrative example we choose a recent application of 
the HPM to a dynamic system of anisotropic elasticity. 

1 Introduction 

In the last years there has been great interest in the application of approx- 
imate variational and perturbation methods that lead to power-series solu- 
tions for linear and nonlinear problems in mathematical physics [1—13]. For 
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example, Chowdhury and Hashim [3] applied the powerful homotopy pertur- 
bation method (HPM) to obtain the Taylor series about x = of the func- 
tions y(x) = e x2 , y(x) = 1 — x 3 /3\, y(x) = sm(x)/x, and y(x) = x 2 + x 8 /72. 
By means of the same method Chowdhury et al [2] derived the Taylor 
series about t — for the solutions of the simplest population models. 
Bataineh et al [1] went a step further and resorted to the even more pow- 
erful homotopy analysis method (HAM) and calculated the Taylor expan- 
sions about x = of the following two-variable functions: y(x,t) = e x2+smt , 
y(x,t) = x 2 + e x2+t , y(x,t) = x 3 + e* 2- *, y(x,t) = t 2 + e xi , and most im- 
pressive: y(x,t) = — 21n(l + tx 2 ) and y(x,t) = e~ tx2 . Zhang et al [12] dared 
to meddle with functions as complex as u(x,t) = —2sech 2 [(x — 2£)/2], and 
u(x,t) = —(15/8) sec h 2 [(x — 5t/2)]. Because of the inherent difficulty in such 
functions the authors pushed the HPM calculation just to first order in t (I 
mean only the linear term of the Taylor series about t — 0). Bataineh et 
al [6] showed great insight and modified HAM to produce MHAM and found 
the Taylor expansions of the functions u(t) = t 2 — t 3 and u(t) = 1 + t 2 /lQ. 
Sadighi and Ganji [11] calculated the Taylor expansions about t — of 
u(x,t) = 1 + cosh(2a;)e~ 41 * and u(x,t) = e 3l ^ x+3t ^ by means of HPM and 
Adomian decomposition method (ADM), and verified that the results were 
exactly the same!!. These authors dared to face the fact that i 2 = — 1. 
By means of HPM Rafiq et al [10] also derived polynomial functions like 
y(x) = x 4 — x 3 , y(x) = x 2 + x 3 and y(x) = x 2 + x 8 /72. Ozis and Agir- 
seven [9] astonished the mathematics and physics community by expanding 
u(x,t) = x 2 e l , u(x,y,t) = y 2 cosh t + x 2 sinh t (and other such functions) 
about t = by means of the HPM. Bataineh et al [5] used HAM to ob- 
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tain expansions about x = for w(x, t) = xe~ l + e~ x , w(x,t) = e x+t+t , 
w(x,t) = e t+x2 and w(x,t) = e t2+x2 . They thus managed to reproduce ear- 
lier ADM results. Although the authors did not state it explicitly in our 
opinion one of their achievements was to prove that the Taylor series have 
exactly the same form in different world locations. 

There are many more articles where the authors apply HPM, HAM and 
ADM and produce results that an undergraduate student would easily ob- 
tain by means of a straightforward Taylor expansion of the model differential 
equations. It is a new trend in mathematical physics and we have just col- 
lected some examples published in only one journal. 

In a recent paper Kogak and Yildinm [13] applied HPM to a 3D Green's 
function for the dynamic system of anisotropic elasticity. The purpose of 
this article is to show the close connection between their results and the 
Taylor series approach. In Sec. [2] we introduce the problem and outline the 
application of the HPM. In Sec. [3] we apply the Taylor-series approach to the 
same problem. Finally in Sec. H]we discuss the results and draw conclusions. 



2 The homotopy perturbation method 

Kogak and Yildinm [13] considered equations of the form 

<9 2 u(x, t) i . 
p ^ ' = Lu(x,t) + f(x,t) (1) 

where u,f 6 R m are vector functions of time and the spatial variables x G 
R n , p G R mxm i s an invertible constant matrix and L a differential operator 
independent of t. 
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They introduced a perturbation parameter p into the equation flTJ 



<9 2 u(x, t) 

} ~dt 2 



P 



Lu(x, t) + f (x, t) 



(2) 



and solved it by means of straightforward perturbation theory (with the fancy 
name of HPM) 

oo 

u(x,f) = £yuW(x,f) (3) 

j=0 

Thus, they obtained 
a 2 u(°)(x,t) 



dt 2 
d 2 u^(x,£) 

dt 2 







Lu^ 1 )(x,t) + f(x,t) ( 5 jl , J = 1,2,... (4) 



Those authors were mainly interested in the Green's function for the differ- 
ential equation (JTJ in which case f(x, t) is a product of delta functions. 



3 Power series approach 

Instead of applying the fashionable HPM we may try the well-known Taylor 
series approach and expand the functions in equation (jTJ) in a Taylor series 
about t — 0: 

oo 

3=0 
oo 

f(x,t) = J>'f,(x) (5) 

3=0 

Thus we obtain a simple expression for the coefficients of the solution: 

u ^+ 2 = ( j + i) 1 ( j + 2 ) p " 1 ( Lu > + fj ) ^ = 0,1,... (6) 
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so that 



t 2 ( *■ \ 
u(x, t) = u (x) + tui(x) + —p~ x (^Luq + f J + 

^r 1 (L Ul + f x ) + . . . (7) 



One easily verifies that all the results shown by Kogak and Yildinm [13] 
are particular cases of equation ([7]) with u = fo = and Lui = — fx so that 
in every case they got the exact solution u(x, t) = tui(x). As an illustrative 
example consider their illustrative example [13] 

Q2 Q2 



L = V 2 = — ,+ 



dx\ dx\ 

f(x,t) = — 2t cos 2 (xi) cos(x 2 ) + 3t sin 2 cos^) 
Mo( x ) = 0, Mi(x) = sin 2 (xi) cos(x2) (8) 

with the exact solution u(x,t) = t sin 2 (a;i) cosf^) that clearly shows what 
we have stated above. 

At this point we want to mention a salient characteristic of the new 
trend mentioned above: the authors fabricate unphysical problems with exact 
solutions and solve them by means of methods like ADM, HPM, and HAM. 



4 Further comments and Conclusions 

The discussion above clearly shows that the paper of Kogak and Yildinm [13] 
is another example of the new trend in mathematical physics that consists 
of obtaining the textbook Taylor series by means of elaborate approximation 
methods. We proved that this standard approach enables one to solve their 
differential equations exactly. However, those authors did not obtain the 
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exact result even at fifth order of HPM. The reason appears to be a wrong 
calculation of the correction of first order. For example, in the introductory 
simple equation they obtained = and = Q(t)tS(x), where Q(t) is 
the Heaviside step function, but this solution does not satisfy the differential 
equation d 2 u^/dt 2 = 5(t)5(x). It is wrong by a factor two and seems to be 
the reason why Kogak and Yildinm [13] did not obtain the exact result at 
second order as predicted by the straightforward Taylor series. However, this 
sloppiness is weightless when contrasted with such magnificent contribution 
to nowadays science and the referees and editors of the journal are content 
with it. 

Another issue is the physical utility of the models considered by Kogak 
and Yildinm [13]. Notice, for example, that the amplitude of the solution 
u(x, t) given by equations (20)-(23) in that paper, which is of the form 
u(x, t) = tui(x) discussed above in Sec. [3l increases unboundedly with time. 
But the new trend will not be deterred by such prosaic considerations. 

The reader may find the discussion of other articles that belong to the 
new trend in mathematical physics elsewhere [14-22]. We recommend the 
most interesting case of the predator-prey model that predicts a negative 
number of rabbits [17]. 

Finally, we mention that a slightly different version of this comment was 
rejected on the basis that "I have determined that it lacks the qualities of 
significant timeliness and novelty that we are seeking in this journal. We 
request you to consider submitting your manuscript in another forum that 
would better suit the material your manuscript covers". After pondering a 
while we realized that we submitted our manuscript a couple of weeks after 
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the appearance of that paper. The new trend is advancing so fast that we 
delayed rather too much. 
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